ARPES spectra and the single particle spectral weight 
of cuprates in the bond-ordered, bond-centered stripe phase 
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The electronic structure and the single-particle spectral density of a stripe array formed by ladder- 
like domain walls (DWs) and by antiferromagnetic (AF) domains of width 2 lattice spacings are 
computed and compared with ARPES spectra from some doped cuprates belonging to the 214 
family of compounds. We assume that bond order is formed on legs in DWs and that the phase 
of the sublattice magnetization changes by tt across each DW. The intensity map plotted in the 
coordination frame momentum-energy reproduces quite well the ARPES spectra obtained at the 
doping level of 15%. We consider this agreement as an argument for a scenario of coexisting bond- 
ordered regions and AF regions in the stripe phase of moderately doped cuprates. 
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INTRODUCTION 



A tendency towards spin and charge ordering in 
cuprates has been observed in results of several neu- 
tron scattering experiments. Low frequency spin fluc- 
tuations observed by many groups at incommensurate 
wave vectors as relatively sharp peeks in the magnetic 
structure factor 0, 0, U have been interpreted as be- 
ing due to stripe fluctuations 0|. Anisotropy observed 
in resistivity [3 also points at stripe formation. Stripe 
order appears to have an impact on phonon mediated 
heat transport Measurements performed by means 
of NMR and NQR techniques demonstrate the emer- 
gence of slow spin fluctuations whose appearance is cor- 
related with pinning of charge- modulations 0- The dis- 
tribution of nearest neighbor bond lengths deduced from 
neutron powder diffraction data ^ Q and measured by 
means of extended x-ray-absorption fine structure spec- 
troscopy 10| agrees with expectations based on a sce- 
nario of lattice response to local charge-stripe order. An 
indication of stripe order may by also found by ana- 
lyzing the shape of ARPES spectra and their evolution 
with doping. For example, it is natural to expect the 
Fermi surface to be flat in the stripe phase in the anti- 
nodal region, near the vectors (±7r,0) and (0,±7r) 4]. In 
addition, quasi-one-dimensionality of the system should 
bring about the depletion of spectral weight in the nodal 
regions, near the Brillouin zone diagonals. Much re- 
search has been done to check if these hypotheses are 
true H m IHIll El m. It seems that some gen- 
eral structure of spectra may be definitely attributed to 
an underlying qua si-one dimensional electronic structure 
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In this paper we will concentrate on doped 
La2-xSrxCu04 (LSCO) compounds, slightly above the 
doping level 1/8 at which, it is believed, a tendency 
towards nanoscale phase separation seems to be ev- 
ident in the the ARPES results El El E3- In 



ARPES intensity maps obtained by integrating over 30 
meV at the Fermi energy the spectra from 15% doped 
La2-x-yNdySrxCu04 (Nd-LSCO) |ll|, patches with 
high intensity may be seen in anti-nodal regions. Appre- 
ciable spectral weight is also detected at the Fermi energy 
in nodal regions. A similar pattern of the spectral density 
at the Fermi energy was predicted by a phenomenolog- 
ical theory of disordered charge stripes and anti-phase 
spin domains (28.] . Unfortunately, this theory does not 
discuss the origin of renormalized hopping terms in the 
effective one-body Hamiltonian and does not explain the 
relation of this renormalization with the underlying mag- 
netic structure. The evolution of the spectral weight as a 
function of doping has been analyzed for the stripe phase 
by means of the cluster perturbation technique (CPT) 
in the framework of the microscopic t-J and Hubbard 
models |H. This theory captures quite well the gen- 
eral trend of this development. Nevertheless, it seems 
that the spectral weight maps derived within this ap- 
proach for the doping level at and above 12.5% do not 
show continuous well developed high-intensity straight 
patches bridging anti-nodal regions. Such structures are 
experimentally seen in nodal regions 13j. In addition, 
that theory does not provide us with much information 
about the magnetic structure of stripes. The remarks 
made above seem to suggest that some understanding of 
the relation between the single particle spectral weight 
of Nd-LSCO at the fiUing level about 1/8 and the for- 
mation of the stripe phase is missing. A phase which we 
may expect to emerge in a natural way in weakly doped 
AFs is a bond-ordered state 29]. Recently, an exact di- 
agonalization of the t-J model (tJM) at a finite cluster 
has been performed to study stripe formation. It has 
been shown that the cluster geometry change, from the 
standard tilted square form of the 20 sites cluster to the 
rectangular one, 5x4 cluster induces the formation of 
a ground state with pronounced stripe-like charge inho- 
mogeneities |30j . The distribution of peaks in the single- 
particle spectral weight, which has been calculated by 
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means of the same method |3l| , resembles experimental 
ARPES spectra from Lai. 28-^Vdo.6'S'?'o. 12(^^04 d- In 
particular, the theory captures quite well the strong dis- 
persion along the (0, 0) - (tt, 0) line. The spreading of 
the quasiparticle-pcaks in the single-particle spectrum 
obtained by means of the exact diagonalization is in 
good agreement with results of an additional calculation 
performed by means of a different method which is the 
bond operator theory for the spin-Peierls phase. Unfor- 
tunately, neither the numerical approach nor the analyt- 
ical method reproduce the flattening of the experimental 
band near the point (tt, 0) and both of them fail to ex- 
plain the emergence of the spectral weight at the nodal 
region, in the form of a straight patch which is observed 
at the Fermi energy in the experiments. Another anal- 
ysis of the spectral weight in doped AFs, based on the 
exact diagonalization of a small cluster concerns the tJM 
with inhomogeneous terms locally breaking the transla- 
tional invariance and the spin-rotational SU(2) symmetry 
[2^ . The pattern of the spectral weight at the Fermi en- 
ergy obtained by means of this method resembles ARPES 
spectra from Nd-LSCO at the doping level 12% and does 
not show enhanced intensity at nodal regions. This en- 
hancement observed in experiments is the manifestation 
of remnant two-dimensionality (2D) in the stripe system. 
Thus, also this theory seems not account for the shape of 
ARPES spectra from Nd-LSCO or LSCO at the doping 
level slightly higher than 1/8. 

Results of inelastic neutron scattering (INS) experi- 
ments indicate that pronounced AF correlations exist in 
the stripe phase On the other hand, the relevance 
of bond-order has been recently demonstrated by means 
of the same experimental method 32]. In a recent work 
[33! we have suggested that the coexistence of bond-order 
with long range AF order may take place in the stripe 
phase. The scenario of nanoscale phase separation is real- 
ized by means of this coexistence. We have shown in the 
framework of this concept that above the doping level 
1/8, when the distance between stripe axes is 4 lattice 
spacing a bond-centered stripe with spin-Peierls order on 
these bonds is more stable than a site-centered stripe. In 
addition it is also known about the site-centered stripe 
that it is more stable than a homogeneous system of holes 
created in the homogeneous AF !£, ^] . The mechanism 
of coexistence between AF order and spin-Peierls order 
in the stripe phase is based on lowering the kinetic en- 
ergy of holes moving without confinement in hole-rich 
stripes which are formed in the form of two-leg ladder- 
like domain walls (DWs) between AF hole-poor domains 
in which the exchange energy decreases 33] . 

In the next section, in the framework of a scenario 
for nanoscale phase separation and coexistence of AF 
long range order with bond-order, we will derive an effec- 
tive tight-binding Hamiltonian describing a quasiparticle 
propagating in such a spin background. Next we will cal- 
culate the part of the spectral function which is accessible 
to measurements in photoemission experiments. Finally 
we will discuss results of the calculation. They seem to 



show characteristic features that may be seen in ARPES 
spectra of Nd-LSCO slightly above the doping level 1/8. 



II. BAND STRUCTURE OF A STRIPE SYSTEM 
WITH COEXISTING BOND AND AF ORDERS 

The tJM in the framework of which we perform the 
calculation is 

H = -T.tA.-^,. + jE(S'S, - 1^). (1) 

The states in which any site is doubly occupied have been 
excluded by definition from the Hilbert space in which 
that model acts. represents pairs of nearest neigh- 

bor (NN) sites. Si and rii denote the operators of elec- 
tron spin and density at site i respectively, t represents 
the hopping matrix element tij between NN sites in the 
square lattice on which the tJM is defined, t'- the hop- 
ping matrix element between second NN sites, and 
t" - the hopping matrix element between third NN sites. 
The rest of tij is zero. We concentrate on the doping 
level about 1/8 at which the distance between axes of 
nearest stripes is 4 lattice spacings and AF correlations 
seem to be of long range 0] ■ Our previous analysis has 
provided convincing arguments that in such a case the 
stripe takes the form of a two-leg ladder-like DW which 
separates hole-poor AF domains 33] . The phase of sub- 
lattice magnetization in domains changes by tt across a 
DW. Each DW is in the spin-Peierls state with singlets 
formed on legs. The underlying spin structure of the 
stripe system at the doping level about 1/8 has been 
presented in Fig^ A natural question arises, what is 
mechanism which gives rise to long range AF correla- 
tions, if AF domains are separated by DWs which consist 
of singlets? The structure depicted in FigQ] emerges in 
the presence of holes only. That presence is induced by 
doping. The creation of a hole gives rise to appearance 
of an uncompensated spin in the DW. This spin may 
be parallel or antiparallel to a nearest spin in the neigh- 
bor domain. The weight of states in which these spins 
are antiparallel, Figl^a), (b), is higher, because such a 
configuration is preferred by AF coupling between NN 
sites. States depicted by Figl^Ja) and (b) are coupled by 
the hopping term in the Hamiltonian. Since their weight 
is higher than the weight of states in which an uncom- 
pensated spin in the DW and the nearest spin in the 
domain are parallel, the hopping term which transforms 
the state depicted by Fig[21(a) into the state depicted by 
FigEIb) ; mediates effective ferromagnetic (FM) coupling 
between sites i and j. Some quantum fluctuations will be 
present in the underlying spin background of the stripe 
system. They may take the form of triplet excitations on 
bonds in DWs and multimagnon excitations in domains. 
FigUin the previous work jSS] contains some examples of 
quantum fluctuations in the spin background. They con- 
tribute a lot to the energy of the system. On the other 
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hand, it seems that quantum fluctuations merely renor- 
malize the shape of quasiparticle dispersion for a given 
spin background and do not change the picture qualita- 
tively. This has been demonstrated in the case of a hole 
propagating in the AF spin background |35l| . a s well as 
in the case of bond ordered two-leg ladders |3(1 Is^ and 
bond-ordered 2D systems "sTl. When it is necessary we 
will take into account the influence of quantum fluctua- 
tions on the distribution of spectral weight. A scenario 
which underlies the calculation which we are going to 
outline next, is based on the assumption that hole mo- 
tion inside ladder-like DWs is governed by the exchange 
of positions between a hole-fermion pair on a bond and a 
singlet on a nearby bond [sBl ■ This exchange is mediated 
by the hopping terms in the initial Hamiltonian Q). We 
also assume that a hole propagates in the AF spin back- 
ground as a spin polaron [Sg ■ In the calculation we take 
into account a simplest form of coupling which moves a 
hole between a bond ordered DW and an AF domain. 
This form of coupling originates in the hopping term of 
the initial Hamiltonian. Before we proceed to construct 
an effective Hamiltonian which describes the motion of a 
quasiparticle in the underlying spin background, we will 
present some formulas which will be useful for that pur- 
pose. A singlet on two sites "L" and "[/", is created in 
the vacuum by the operator 



belongs to a bond occupied by a singlet and vice versa 



(2) 



a° is the two-dimensional identity matrix and cr", a = 
X, y, z, are Pauli matrices. The summation over repeating 
Greek indices is assumed. Three operators which are 
components of the vector t\^^ create three triplet states 
on sites L and U, 



"LU 



V2 



(3) 



A formula which wc will often use later is 
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it[,y,cr„^c[„)|0). 



(4) 



The left side of Q together with the right-side term 
which contains a singlet at the pair of sites L' and U' 
represent the exchange of a singlet and a hole-fermion 
pair between two bonds, which is mediated by the hop- 
ping term in the initial Hamiltonian. As we have already 
stated before, we will neglect the creation of triplets on 
bonds, however this process also takes place during hole 
hopping. The creation of a triplet is represented by the 
second term at the right side of Q . We also need formu- 
las which represent the hopping of a hole from a single 
site i which belongs to an AF domain to a site which 
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When we start to analyze contributions from quantum 
fluctuations in the ground state of the system to the spec- 
tral function we will need foUowings formulas 



1 
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t J 



)|o). 



(7) 
(8) 



We begin the derivation of the effective Hamiltonian de- 
scribing hole propagation in the spin background de- 
picted by Fig^ with outlining the mechanism of hole 
propagation inside an AF domain. A hole created and 
moving in the Neel background shifts spins between dif- 
ferent sublattices and creates defects in the AF structure, 
FigEl (a')-(c). Such a process gives rise to an increase of 
the Ising energy. This rise is roughly speaking propor- 
tional to the length of a path along which the hole has 
traveled, which means that a tendency towards hole con- 
finement appears |^ . In order to take into account such 
a tendency we will analyze hole motion in the framework 
of a basis which consists of states representing holes con- 
fined in the AF background by linear defects (strings) left 
behind by moving holes on their way. We call these states 
spin polarons. A wave function representing a confined 
spin polaron at a site i in an AF domain is a combination 
of states which are created from the state Ci_m)|iV) by 
the NN hopping term when a hole created at the site i 
starts to move, 



(9) 



\N) is the Neel state in the domain, iVi) denotes a state 
obtained by hopping along a path Vi of a hole created at 
the site i. ajjp.) is the amplitude of this state. We have 
assumed for simplicity that cticPi) depends solely on the 
length /(T'i) of the path Vi- The length of a path or of 
a string state is defined as the number of hops needed 
to form a given string state from a state representing a 
hole created in the AF spin medium. At the first stage 
of the analysis we take into account only the hopping be- 
tween NN sites because t 3> J,t',t". Processes related 
to hopping between further neighbors and processes re- 
lated to swapping antiparallel spins by the transversal 
term in the exchange interaction will be considered later 
as a perturbation. A hole moving inside a domain may 
make its first step in {z — 1) directions. 2=4 is the co- 
ordination number of the square lattice. There are in 
principle (z — 2) direction choices of each next hop, if the 
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hole moves without retracing inside the domain. On the 
square lattice there are z choices for the direction of the 
first step and (z — 1) for the direction of further hops 
during the non-retractable motion. These number gets 
reduced by one for the domain formed by two chains of 
sites. Thus if we neglect some details, as for example 
path crossing, we may write, 

(M'.l*.) = + - 1) ^(^ - ly-'al. (10) 

Each prcfactor in IjlOl) at the square represents the 
number of different paths with the length /i. We calcu- 
late the energy of the spin polaron state ei, which 
is given by the expectation value of a trial Hamiltonian 
i/o, (^'il-ffol^i): with the assumption, that the motion 
of a hole which has started from the site i, is restricted 
to the interior of the domain, and that the contribution 
from the interaction term is restricted to the Ising part 

+ 2(z-l)^(z-2)'^a^a^+ii. (11) 

The first term in Hll|l basically counts the number of 
"broken bonds" , that are not occupied by a pair of an- 
tiparallel spins, in which case the Ising contribution to 
the energy of that bond is higher by J/2 than in the 
case if it were occupied by a pair of antiparallel spins. 
The second term in (|11|) is the contribution from the 
hopping operator to the spin polaron energy. The pref- 
actors appearing in this term represent the number of 
paths with a given length multiplied by the number of 
directions in which these paths may be extended. The 
appearance of the factor 2 is related to the fact that the 
hopping which couples paths of length [i and ^ -I- 1 may 
take place forth and back. The values of these paths 
can be found by minimizing under the con- 

straint i^i\R\^i] = 1. After we have constructed spin 
polarons which are formed in AF domains we are able to 
present the full basis of single-particle states. The under- 
lying spin background which has been presented in Fig^ 
plays the role of the vacuum \Q) for hole-like quasipar- 
ticles which propagate in this background. has been 
obtained by acting on the absolute vacuum |0), in which 
no particles are present, with a product of operators like 
s]^u creating singlets on bonds connecting sites L and 

U and operators creating spins in domains accord- 
ing to the pattern shown in Fig^ New fermionic pos- 
itively charged operators create single particle hole- 
like states from the vacuum for holes The action of 
the operator /lio-l^t/a-) '^^ 1^)' where the site L{U) be- 
longs to a ladder- like DW, exchanges the operator s]^u in 
the product defining by the operator clj^{c\^^), which 
means that instead of a singlet on the bond connecting 



sites L and U there is a hole on the site L{U) and spin a 
on the site U{L). The action of the operator h]^ on a site 
i which belongs to an AF domain, creates a spin polaron 
I^E'i) in that domain. The spin polaron is a combina- 
tion of some states, the amplitudes of which are given by 
prefactors a^. These states include Cia\^) and states ob- 
tained by applying consecutively the NN inside-domain 
hopping term to the state Ci^\Q). During this process, 
the hole hops without retraces and fi is l{Vi), the length 
of the path Vi, along which the hole has traveled to form 
a given component state of A label which we will 

use to mark the fermionic operator creating either a bond 
hole or a hole-like spin polaron at a give site is {"^'^) ■ rn 
refers to the column number of the unitary cell to which 
belongs the site where the hole-like particle has been cre- 
ated, n refers to the row number of the unitary cell, and 
i (column), j (row) are indices representing the position 
of that site inside the unitary cell, i and j run from to 
7 and from to 1, respectively. Now we start to explain 
with some details the origin of contributions to the effec- 
tive Hamiltonian for a single quasiparticle propagating in 
the spin background depicted by Fig ^ which is an exem- 
plification of coexistence between AF and bond orders. 
We concentrate on the case of a spin-up quasiparticle. 
Within the approximation which we use, operators 
create in the vacuum |r2) eigenstates of the unperturbed 
Hamiltonian iJg, 

i/o = E s^s, - 1 y: + J E is!s^ - ^) 

(12) 

where are pairs of sites on which singlets depicted in 
Fig^have been formed and are pairs of NN sites 

belonging to AF domains. Action of Hq is restricted to 
the space containing states in which none of the sites is 
doubly occupied. This Hamiltonian contains the Ising 
part of the exchange energy of links inside domains. It 
also drives hole hopping between NN sites inside each do- 
main. Since all matrix elements which may give rise to 
deconfinement of a hole, have been by definition removed 
from Hq, polaron states l^f^) are its eigenstates. Pro- 
cesses which bring about deconfinement of holes will be 
treated at the latter stage of the calculation as a pertur- 
bation. The exchange energy and the hopping between 
sites belonging to bonds occupied by singlets in Fig^are 
the only contribution to the Hamiltonian Hq from DWs. 
Terms in the Hamiltonian of the tJM which couple a site 
belonging to a ladder-like DW with a site belonging to 
a domain do not contribute to Hq. It is clear that the 
vacuum state and the single particle states hl^\n) are 
eigenstates of Hq. Within the first order approximation 
the on-site energy of a quasiparticle with spin up created 
at the site ('JJ''") which belongs to a DW is 2J. From 
now on, the reference value of the energy is the energy of 
the vacuum state The contribution from a destroyed 
singlet to the on-site energy of a hole-like quasiparticle 
created at the site (™'") is (3/4)J. A hole-like quasipar- 
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tide occupying the site ("J^'f ) with spin up is by definition 
the same as a single fermion with spin up which occupies 
the site ("^'i^) belonging to the bond ('o'(^)-(™'") • Since 
spins on the sites (™g') and (™f q ") are in this case paral- 
lel, the lowest order contribution to the exchange energy 
of the bond between these two sites additionally increases 
by J/4 compared to the contribution from this bond in 
the vacuum state In the presence of a hole, the con- 
tribution from the potential term — j) '^i'Tij in the 
initial Hamiltonian of the i JM is higher by J. By adding 
all partial contribution we get the value 2J of the to- 
tal on-site energy of the spin-up quasiparticle at the site 
(™'"). The same on-site energy have spin up quasipar- 
ticles created in DWs at all sites which are NN of sites 
in domains occupied by spins pointing down. Finally we 
deduce that in the effective single-particle Hamiltonian 
appears a term 

+ ^b)T^(lo')T+^b)T''(lo")T]- (13) 

The fermionic operators ^|m,,i-j and h^m,n^ transform the 

underlying vacuum into the single particle state de- 
scribed above and the single particle state into \ft), re- 
spectively. The notation for indices which we use seems 
not to be very short, but such a form of it is basically 
unavoidable because the elementary cell of the underly- 
ing spin background, Fig^ is rather big. That notation 
also helps to trace easily at the map which is Fig^ the 
results of hopping events mediated by the Hamiltonian 
Hgf f . An analogous term will appear in the Hamilto- 
nian representing a propagating spin-down quasiparticle, 
but with a different set of indices labelling sites in the 
elementary cell. The creation of the spin down quasipar- 
ticle at these sites, (0,0), (1,0), (4,1), (5,1) induces the 
formation at the ladder leg of an uncompensated spin, 
the direction of which is parallel to the direction of the 
nearest spin in the domain. For the sake of simplicity we 
will concentrate in this paper on the propagation of the 
spin-up quasiparticle. Such a simplification is possible 
because despite the breakdown of the time reversal sym- 
metry by the underlying spin structure, the energy of the 
spin-up and the spin-down quasiparticles is degenerate. 
Later we will user this observation in the calculation. The 
on-site energy of the hole-like quasiparticle is lower for a 
group of sites by J/2 then for sites to which the contribu- 
tion (|13() refers because spin of the bond-fermion which 
appears after a hole has been created at a bond initially 
occupied by a singlet may be antiparallel to the nearest 
spin in one of domains, 

S2H.ff = 2^E['^(-„")T^"o")T+^(7o")T^(to")T 

m.n 

+ ^K;r)T^Tr)T + ^(lT)T^'?:r)T]- (i^) 



Let us concentrate now on the on-site energy of quasipar- 
ticles created at sites belonging to domains. A spin-up 
hole-like spin polaron can by definition be created ex- 
clusively at sites which have been initially occupied by 
a spin-down fermion. An obvious contribution to the 
on-site energy of a hole-like quasiparticle energy in do- 
mains is El, the minimum value of the matrix element 
obtained under the constraint (^'^j^'i) = 1. During 
the construction of spin polarons we have considered only 
hopping between NN sites, which is governed by the hop- 
ping with the highest integral t. We have also neglected 
some "high order" processes related to path crossing or 
to the action of the XY term in the Heisenberg model. 
In the analysis which we start now, we will discuss in 
the framework of the first order perturbation theory the 
contribution of some neglected processes to the effective 
Hamiltonian iJe//- States depicted by Figs. 4(b) and (c) 
are string states of length 1 and are components of the 
spin-polaron wave- function j^'^) at the site i. A single 
vertical or horizontal hop of the hole created in the AF 
background at the site i gives rise to the states depicted 
by Figs^l^b) and (c) respectively. The hopping term to 
next nearest neighbors (NNN) which we treat as a per- 
turbation couples states represented by Figs. 4(b) and (c), 
which brings about a contribution to the matrix element 
(^'ilJfil^'j), where Hi = H - Hq. This contribution is 

7i=4t'a2. (^5) 

We recognize in (|15|) a product of amplitudes with which 
string states of length 1 appear in the definition © of the 
spin polaron state. The factor 4 which appears in H15() 
is related to the fact that hopping which couples states 
(b) and (c) may occur in both directions and that analo- 
gous coupling takes place between the state depicted by 
Fig^c) and the state depicted by Fig^d). The state 
depicted by Fig^Jd) has been obtained by a single down- 
ward hop of a hole created in the AF domain at the site 
i. The contribution 

72 = 2t"al (16) 

to {^i\Hi\'^i) originates in a similar way with coupling 
between states depicted in Fig^c) and (d) by the hop- 
ping term to third nearest neighbor (TNN) sites. The 
coupling between longer string states which are compo- 
nents of the same spin-polaron wave-function may also 
contribute to the renormalization of the spin-polaron on- 
site energy. To be more specific, the coupling between 
states as depicted by FigsQJe) and (f) and between their 
refiections in the horizontal line running through the site 
i gives rise to the correction 

73 = 2t'[2al + {z-l) ^(z - 2r-'al] (17) 

to {'i'i\Hi\'i'i), while the coupling between the state de- 
picted by Fig0Jg) and its reflection in the horizontal line 
running through the site i together with some similar 
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process in which longer strings are involved brings about 
the correction 



formed in the underlying spin background 



74 = 2t"[a2^(z- 



fi—3 



(18) 



The states represented by Figs^e), (f) and (g) have been 
obtained by means of three different sequences of hole 
moves. The hole has been created at the site i and these 
sequences are upwards- upwards, upwards-left and left- 
upwards respectively. Longer strings pinned to the site 
i may be coupled in a very similar way, which also gives 
rise to a change of the on-site energy. This change has 
actually already been incorporated into parameters 73 
and 74. Finally, by collecting all terms we may infer that 
in the effective Hamiltonian -ffe// appears the following 
term referring to the on-site energy of spin-polarons 

SsHeff = (si + 71 + 72 + 73 + 74 + J/4) 



(-)T^To")T-^"(7;r)T"C?:DTJ 



(19) 



The additional term J/4 in the prefactor is related to the 
fact that the contribution from the contact interaction - 
JniUj /A between a site z in a domain and a site j in a DW 
is higher by J/4, when a hole is created in the domain 
at the site i. This interaction was neglected, when we 
were calculating the eigenenergy Si of the spin polaron 
and needs to be taken into account now. (|19|) is the last 
on-site term in the effective Hamiltonian. 

The NN hopping integral for the quasiparticle, the 
propagation of which we describe in this paper is the 
same as the bare hopping integral provided that the 
sites between which the quasiparticle moves belong to 
the same initially ordered bond inside a DW 33]. Thus, 
the term which describes the quasiparticle hopping be- 
tween NN sites, on which singlets have been formed in 
the underlying spin background, takes the following form 

rn.n ' 



(20) 



The NN hopping integral for a pair of sites, which belong 
to different singlets in the underlying insulating state is 
t/2 instead oi-t. This change of sign and this reduction of 
size may be deduced from the form of the first term on the 
right side of Eq.Q. Therefore, we may write the follow- 
ing expression for the contribution to the effective Hamil- 
tonian describing quasiparticle hopping between NN sites 
belonging to different bonds on which singlets has been 



65H, 



t 



(-)T^-)T+^— ^)T^-)T 



+ '^(-)T^(":o")T + ^r.':r^)T'^C::.")T 



'^(-)T^Tr)T + "("xr^)T"(i'^)T 



(21) 



The value of the NN hopping integral, between a site 
which belongs to a domain and a site which belongs a 
DW, can be inferred from the first part of the right side 
in Eq.@, 



(IT)T (i,i)T ("'oOT (4:0)1 



]+H.c.}. 



(22) 



Since the time reversal symmetry and the translational 
symmetry are broken inside AF domains, the propagat- 
ing quasiparticle can not move between different sublat- 
tices and any term related to NN hopping inside domains 
is not generated in the effective Hamiltonian iJe// ■ NNN 
quasiparticle hopping inside DWs and hopping between 
sites which belong to a DW and a domain is a first order 
process that is mediated by the NNN hopping term in 
the bare Hamiltonian. The explicit forms of related con- 
tributions to Heff may be deduced from formulas l@J, 
© and 



SrHeff 



t' 



(23) 



+ '^riT^)T^To")T + ^(7,")T^":r)T 

- -^"oE{['^;-"--)t'^(-„")t 

+ '^rTr^)T^";o")T + ^b)T^":")T 



The task of finding the formula for the term describ- 
ing NNN hopping inside AF domains is little bit more 
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tedious. For example the coupling by the XY term in 
the Heisenberg model between string states depicted by 
Fig^f ) and (h) gives rise to hopping terms in H^f / which 
shift a spin polaron from the site i to the site j and 
vice versa. We analyze now the XY term, because it has 
been neglected during the first stage of the analysis, when 
quasi-confined spin polaron states have been constructed. 
Since the string state depicted by Fig0Jf ) is a component 
of the wave function j^I^i) for the spin polaron created at 
the site i and the string state depicted by Fig^fh) is ev- 
idently a component of the wave function for the spin 
polaron created at the site j, the coupling between these 
components gives rise to the coupling between spin po- 
laron states IvPi) and This brings about a contribu- 
tion to the matrix element liyl^Pi) and to the hopping 
term in the effective Hamiltonian -ffe//- It is probably 
useful to remind now that H^ff is an approximation to 
the initial Hamiltonian H in This approximation 
is expressed in terms of the operators creating and an- 
nihilating spin polarons. Also the coupling between the 
states depicted by Fig^g) and (h) contributes to the 
hopping term between the sites i and j in the effective 
Hamiltonian H^ff- This coupling is mediated by the XY 
term. We have already discussed the action of the XY 
term which by removing two defects in the AF structure 
transforms string states of length 2, tails of which are 
pinned to the site i, into the state representing a hole 
created at the site j in the AF ordered domain. By a 
string tail we mean its end opposite the end at which sits 
a hole. The XY term may also transform a state repre- 
senting a hole created in the domain at the site i into a 
string state of length 2 pinned at the tail to the site j. 
This additional coupling between components of the spin 
polaron states at the sites i and j doubles the value of 
the hopping integral between these sites in the effective 
Hamiltonian. Thus, after a little thought we may deduce 
that the integral for the NNN hopping inside domains is 



2J^(z 

11=2 



2y 



(25) 



The first term in the sum presented above refers to cou- 
pling between strings of length and 2. We have just 
outlined its origin in detail. Other terms appear in the 
sum (|25|l because longer strings, which are created when 
a hole moves further from the site j in FigsQ^f), (g) and 
(h), are also coupled by the XY term in the Heisenberg 
model. The hop left of the hole depicted by Fig^Jg) 
gives rise to the string state of length three, Fig0fi), 
which is pinned to the site j. Since we neglected the 
possibility of path crossing when we were constructing 
the quasiconfined spin-polaron states, some corrections 
need to be made now. For example we did not consider 
before, that by applying the NN hopping term to the 
state depicted by Fig^i) we may create the state de- 
picted by Fig^j). Since the latter state is a string-hke 
component of the spin-polaron at the site j obtained by 
hopping downward and left of a hole created at that site, 
we deduce that the process described above generates the 



NNN hopping term in the effective Hamiltonian with the 
amplitude 



T2 — 2ta3a2- 



(26) 



The factor 2 originates with the fact that the motion of 
a hole around a plaque in the square lattice may take 
place clockwise and anti-clockwise. We also recognize 
and a2 as amplitudes of strings, which have the length 
3 and 2, respectively. Exactly such strings which are 
components of the spin-polaron states at the sites i and 
j in Figs^i), (j) are coupled by the NN hopping term 
in the initial Hamiltonian. The NNN term in the ini- 
tial Hamiltonian generates coupling between states rep- 
resenting "bare" holes created in the AF background of 
domains. Since these states are also string components 
with the length of some spin polaron wave functions, 
coupling between the latter is also generated. An exam- 
ple of such coupled string states are Figs0Jh) and (k). 
The contribution to the NNN hopping integral in H^ff 
is 



T3 = t'al- 



(27) 



The NNN term in the bare Hamiltonian also couples 
states depicted by FigsQff) and (1) which are string-like 
components with length 2 of spin-polaron states at the 
sites i and j, respectively. The coupling amplitude is 

r4 = 2t'al (28) 

The appearance of the factor 2 is related to the fact that 
states which are created when holes move between the 
sites i and j in opposite directions around the plaque, 
than in the case of the states depicted by Figs^Jf) and (1) 
are also coupled by the bare NNN hopping. The state in 
Fig2Jl) has been obtained by hopping left and downward 
of a hole which has been initially created at the site j. 
By collecting all contributions, which we have discussed 
above, we see that the new term in to i?e// is, 

SgHeff = in +T2+T3+ T4) 



(-)T^IT)T^"(-7^)T"(":^)T 



H.c.}. 

(29) 

By means of a similar analysis, as for the hopping be- 
tween NNN sites we may find the TNN hopping term in 
the effective Hamiltonian. For the operator representing 
the quasiparticle hopping inside DWs we get 



(-:o+^)t"("o")t 
r)T 



''ri:o+')T''(To")T 



r)T 



C".':o+^)t"(Io")t 

Ot+'^('"-+i)T^Io")T 
+ H.c.}. (30) 
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The TNN hopping of the quasiparticle between domains 
and DWs is governed by the following term 

t" 

+ ^(-;r)T^(o;r)T + ^(-)T'^(Io)T 

+ /.}„_,„)^/^(..),+/.;_)^/^(.^ 

+ hy^fy^^..)^]+H.c.}. (31) 

To TNN hopping term between sites in domains con- 
tribute: a) the mechanism which is based on the short- 
ening of strings by the action of the XY term, b) the cou- 
pling of strings with length by the TNN hopping term 
in the initial Hamiltonian and c) the exchange between 
the head and the tail of a straight string with length 2. 
The last process is mediated by the TNN hopping term 
in the initial Hamiltonian 

The contributions to the hopping integral in these 
three cases are 



M=2 



r7 = t"al 



(32) 

(33) 
(34) 



Since the origin of these couplings is the same as for the 
NNN hopping term in -ffe// we do not discuss them in 
detail. Thus, the TNN hopping operator, which is the 
last contribution to i?e// discussed by us at the approx- 
imation level, that we have assumed, takes the form. 



('"c+^)T"(Tr)T 



^rTo+^)T^Io")T+'^('-+^)T^'S,o")T 



\"-^)t"(?;Dt 



]+H.c.}. 



(35) 



Figl^l depicts the electronic structure which we have ob- 
tained by solving the Hamiltonian Hej f . It represents the 
energy dispersion E(p) of all bands along the line (0, 0)- 
(tt, 0) -(tt, 7r)-(0, 0) and the line obtained by performing 
the rotation by 7r/2 around the point (0, 0). Such a com- 
bination of dispersion curves is justified by conditions 
in which experiments are performed, because it seems 
that stripes in some regions of the sample may run ver- 
tically, while in other regions they may run horizontally. 
The reason for this mixing can be, for example, sam- 
ple twinning. Nondispersing parts of bands are related 
to some obstacles for hole propagation in the directions 
perpendicular to stripes. That kind of motion seems to 
be blocked, which is clear because i?e// by definition can 
not mediate quasiparticle motion occurring exclusively 
in the direction perpendicular to stripes. It seems that 



the lack of the energy dispersion in some directions ac- 
counts for straight patches of the spectral weight which 
appear in the density maps obtained by means of ARPES 
measurements '40^. The coordination frame applied to 
draw these maps is the same as we have used to obtain 
FigEl (momentum-energy) , while the second derivative of 
ARPES spectra plays the role of a density-like parameter 
which marks regions with the high spectral weight ^27j. 
We shall make a more detailed comparison with ARPES 
spectra, after we calculate the spectral weight by means 
of the approach applied by us. That approach is based, 
as we have already mentioned several times, on a combi- 
nation of bond and string formalisms. 



III. SINGLE PARTICLE SPECTRAL 
FUNCTION IN THE STRIPE PHASE WITH 
COEXISTING BOND AND AF ORDERS 



ARPES probes the one-particle spectral function 
A~{'k,uj). We neglect in the further analysis of ARPES 
spectra the influence of temperature. We also assume 
that the description of the photoelectric effect in terms 
of Fermi's golden rule is sufficient and finally we omit, in 
the calculation, electromagnetic dipole matrix elements 
between the wave function of photoelectron and the wave 
functions of the electrons in the initial states. Further- 
more, we apply the single particle approximation, in the 
calculation of the spectral function A~(k, w), which is 
natural, because the Hamiltonian -ffe// by definition does 
not contain the interaction terms. Since electrons are 
emitted by the photoelectric effect, at T = information 
is gathered only about the one-electron removal part of 
the spectral function, which takes the form 



Af-l| 



(36) 

Due to the formation of the stripe structure the shape of 
which is determined by the underlying spin background 
depicted by FigQ] the IBZ gets reduced by a factor of 
8 in the horizontal direction and by a factor of 2 in the 
vertical direction. The Hamiltonian H^f f may be written 



in the diagonal form in terms of operators h 



t 

k,T,a- 



and 



their Hermitian conjugates, which are determined by the 
form of Hamiltonian eigenstates, 



h 



t 



1 



(»j),T,' 



("^)T 



(37) 

where k/j belongs to the reduced IBZ, a labels the band 
number, d = 8 is the length of the elementary supercell, 
if = 2 is the width of the elementary supercell and Nd x 
Lw is the system size. Thus, within the single-particle 
approach, the one-electron removal part of the spectral 
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function is approximated as 
A-{'k,uj) ^ ^(5(t^ + ek,a) 



I X! X! X! K.Tnod-KR,{i.j);ha^ 
n.l i,j i' ,j' 



ik(ncZ,/i(j) 



X e 



(r!|/i(„,,)^C(^o,o^)jr!)|2, (38) 



where kmodK^ denotes the vector k reduced to the IBZ 
of the superlattice, the elementary ceh of which is de- 
picted by Fig^ and £k,Q is the energy of the a-th band 
in FigEl In order to evaluate we need to find matrix 
elements, 

Af^,i.i-jj^ 0,0 -J = (57|/i^n,!'j^c^ 0,0^-j J57). (39) 

A scheme showing how to do this in the framework of 
the spin polaron (string) approach was developed before 
|4l|. For example, the removal of the spin down electron 
at the site (2,1) in the elementary cell of the underlying 
spin background depicted by Fig^ gives rise to a state 
which is component of the wave function for the polaron 
created at this site. Thus the matrix element H39|l for 

= ihj) = (2,1) and {n,l) = is ao- In the ap- 
pendix we list all sites labelled by the pairs of numbers 

{n,l), for which a nonvanishing matrix el- 

ement (I39II exists. The removal of a spin down electron 
from a site which belongs to ladder-like DWs in the state 
depicted by Fig^ gives rise to the hole-like quasiparticle 
h at this site. Since this removal may be mediated by 
the electron annihilation operator c in H39(l , a nonvanish- 
ing diagonal matrix element M = — 1/\/2 is generated. 
Quantum spin fluctuation which are generated in the un- 
derlying spin background schematically depicted by Fig^ 
can not be neglected when we evaluate the spectral func- 
tion. Within the first order perturbation theory the ad- 
mixture 5\4>) of quantum fluctuations to the groundstate 
\4)o) of Hq brought about by the perturbation Hi is 



En — En 



(40) 



where l-^n) are excited eigenstates of Ho with energies 
En and Eq is the groundstate energy. The action of the 
exchange interaction between sites depicted by Fig|Sfa), 
L{U) belonging to a DW and belonging to an AF do- 
main, may transform the singlet on sites L and U into a 
triplet. Within our approach we treat this part of the ex- 
change interaction as the perturbation. The transformed 
state is an excited eigenstate of Hq with the energy J 
and also a quantum fluctuation that, according to the 
formula (|40|l , contributes to the underlying vacuum state 
|f2), about which we assume that it is an eigenstate of the 
Hamiltonian QJ. jfi) is schematically depicted by FigQ] 
in which quantum corrections have not been taken into 
account even within the framework of the first order per- 
turbation theory. The electron annihilation operator c 
that acts on a site belonging to a bond, which has been 



excited to the triplet state, transforms it into the state 
representing the hole-like quasiparticle h created at that 
site. This process gives rise an addendum to the ma- 
trix element (|39|) . This addendum may be attributed 
to the existence of quantum corrections to the vacuum 
state For example the value of this addendum is 

-l/(2\/2) for the matrix element (|39|l labelled by the in- 
dices (i',/) = = {n,l) = (0,0). The XY part of the 
exchange interaction which we treat as a perturbation, 
creates in an excitation that takes the form of two 
NN spins turned upside down, Fig|SJb), with respect to 
the underlying spin structure, Fig^ According to the 
recipe this excitation contributes a correction to the 
underlying vacuum state If the annihilation oper- 
ator c removes the spin down fermion at the site j in 
FigEIb), the configuration depicted by FigEl^c) will be 
created, which is a string state, a component of the wave 
function [ipi) for the spin polaron at the site i. Since 
this spin polaron state is by definition created by the 
fermionic operator a contribution -ai/4 to some di- 
agonal matrix elements of the form (|39|l is generated. 
We discuss now the last category of processes which give 
rise to new terms in the sum that appears in H38|l . The 
electron removal from the site k, Fig^fb), gives rise to 
a string state of length 2, Fig^d), which is a compo- 
nent of the wave function for the spin polaron at the 
site i. Since the spin-polaron state is created by the op- 
erator , a nonvanishing matrix element is generated. 
Its value is -a2/4, which gets multiplied by a factor of 
2, because there are two strings of length two, connect- 
ing sites i and k. By collecting all contributions to the 
sum which appears in H38|) the one-electron removal part 
of the spectral function j4^(k, w) can be evaluated. We 
will apply in the numerical evaluation some Lorentzian 
broadening of the Dirac delta function, which is justified 
because experimental measurements, with which we are 
going to compare our results, have finite resolution and 
also some averaging procedure is often applied to present 
experimental data. 



IV. NUMERICAL EVALUATION OF THE 
SPECTRAL WEIGHT, COMPARISON WITH 
RESULTS OF ARPES EXPERIMENTS AND 
CONCLUDING REMARKS 

It seems that a most complete set of data to compare 
with our theoretical analysis, provide ARPES measure- 
ments of the LSCO and Nd-LSCO systems not for the 
doping level x = 0.125 but for the doping level x = 0.15 
\l2L fla . l40l |. We will present our results in a way 
similar to the method of presentation, which was used 
in experimental papers reporting these measurements. 
The NN hopping parameter t defines a unit in which 
energy is measured. We choose J/t = 0.4, t' /t = —0.1, 
t" /t = 0.05. Similar parameters have been applied in 
a recent theoretical analysis of ARPES spectra in the 
2D tJM, based on the exact diagonalization of the 5x4 
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cluster and in a separate calculation performed by means 
of the bond operator formalism applied to the columnar 
bond order underlying spin structure j^. This coinci- 
dence helps to make comparison between results of the 
calculation presented in our paper and results of the ear- 
lier analysis. The choice of the Hamiltonian parameters 
is ever to some extent arbitrary. On the other hand a 
theoretical analysis of ARPES spectra from LSCO sys- 
tems, based on the tight binding approach indicates that 
the ratios \t' /t\ and \t" /t\ are lower for these systems 
than for other members of cuprate family |42j| . Parame- 
ters suggested by the authors of the analysis based on a 
tight binding approach 3.re basically the same as pa- 
rameters applied in our calculation. The position of the 
Fermi energy in the band structure at the doping level 
15% has been determined in our calculation by counting 
the number of hole-like states. It has been marked as a 
narrow straight line in Fig|H| Fig[7| depicts the intensity 
of the one-electron removal spectrum function presented 
in the coordination frame momentum-energy. Contri- 
butions from vertical and horizontal stripes have been 
summed in order to account for presumed coexistence 
of these structures in different parts of the sample. In 
some agreement with the experimental result obtained 
at the doping level 15% |23| we notice in the calculated 
spectrum a strongly dispersing band between points (0, 0) 
and (tt, 0)/(0, it) which approaches the Fermi level, uj — Q, 
near the zone boundary, where it joins a flat patch formed 
by the region of high spectral intensity. After passing the 
anti-nodal region as a straight narrow strip, the band-like 
region of high intensity disappears somewhere between 
(tt, 0)/(0, tt) and (7r,7r) points. In the calculated spec- 
trum we also notice two cusps near k — (7r/2,0) and 
k = (7r,7r/2) which are absent in the experimental spec- 
trum. This discrepancy may be attributed to the fact 
that to this high-intensity patch which looks like a single 
band in the plot, actually contribute two of many bands 
which may be seen in Fig|3 The band which appears 
in the results of the exact diagonalization performed for 
the i JM and the results of the calculation based on the 
scenario of the columnar spin-Peierls order jSlJ has also 
the strong dispersion but does not flatten in the anti- 
nodal region. Despite that neither the spectral weight of 
the quasiparticle propagating in the bond ordered spin 
background nor the spectral weight of the quasiparticle 
propagating in the AF background 43] resembles ARPES 
spectra from doped LSCO and doped Nd-LSCO in the 
stripe phase, spectral properties of the model in which 
these two phases coexist gives rise to qualitative agree- 
ment with experimental data. This agreement may be 
attributed to the fact that the band structure of the 
model of nano-scale phase separation which we discuss 
here has some features which are present in band struc- 
tures of both pure homogeneous phases. For example in 
our results we actually see a remnant of a strip formed 
by high intensity which takes the shape of a strongly 
dispersing band along the line between points (0, 0) and 
(tt, 0)/(0,7r) and resembles the band that is formed in 



the columnar bond-ordered phase. This band-like struc- 
ture actually does not flatten and reaches the maximum 
near the points (0, 7r/2)/(7r/2, 0). On the other hand it 
is likely that the presence of this maximum may be at- 
tributed to the brute force method of sewing the bond- 
ordered parts of the system with the parts which are 
AF-ally ordered. One thing is for sure, that the lack 
of the straight high intensity patch at the anti-nodal re- 
gion in results of the theory based on the scenario of 
the columnar spin-Peierls phase demonstrates that the 
spectral function calculated within this picture does not 
agree in some details with measured ARPES spectra and 
suggests that it is necessary to take into account also the 
long range AF correlations in order to formulate a theory 
which accounts for the spectral properties of the cuprates 
in the stripe phase. In the region between the zone center 
and antinodal points (tt, 0)/(0, tt) similar agreement has 
been observed between the experimental data and results 
of a calculation based on the dynamical mean field theory 
(DMFT) applied to the Hubbard model 

Along the line connecting the points (0, 0) and (tt, tt) 
both in ARPES spectra and in our results we see a band 
in the vicinity of the Fermi surface. This band has a 
maximum at the points (7r/2,7r/2), bends downward and 
disappears near the zone corner. A similar feature may 
be observed in the results of an exact diagonalization 
performed for the t-t'-t"-J model and in the results of the 
calculation performed within the scenario of the bond- 
ordered columnar phase The DMFT of the stripe 
phase in the Hubbard model at the doping level 15% gives 
rise to a slightly different result, namely it seems that the 
band crosses the Fermi level near the point (7r/2, 7r/2). 

The intensity map of the one-electron removal spec- 
tral function (k, oj) at the Fermi energy obtained in 
the framework of the scenario which assumes coexistence 
of bond and AF orders is depicted by Fig|S| FiglSfa) 
refers to results of the calculation in which we have 
assumed that bond order is formed on legs in ladder- 
like DWs. That figure also shows some agreement with 
ARPES data from LSCO and Nd-LSCO O obtained 
for the doping level 15%. We see well developed spectral 
weight in the nodal regions. These regions are bridged by 
high-intensity continuous almost straight high-intensity 
patches. It seems that the agreement between the exper- 
imental results and the results of the theory based on the 
mixture of the bond formalism and the spin polaron ap- 
proach is better in this respect than the agreement with 
results of previous calculations based on the CPT ap- 
plied to microscopic models which are the t JM and the 
Hubbard model 20, 25]. 

FiglHIb) depicts the spectral weight at the Fermi en- 
ergy obtained in a separate calculation for the under- 
lying spin background with bond order on rungs. We 
know from the results of a previous paper [s^ that such 
a structure is less stable and do not expect much simi- 
larity with experimental results. Such a lack of similar- 
ity may be seen, indeed. For example the high-intensity 
patches does not form a shape resembling the Fermi sur- 
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face obtained by means of calculations based on the local 
density approach. Such a shape may be observed both 
in experimental spectra and in FiglHIa). In FiglHl apart 
from patches of high spectral density we notice additional 
regular structures formed by regions of enhanced inten- 
sity. It seems that the origin of thousc structures may 
be attributed to simplicity of our approach within which 
fluctuations of the underlying spin structure depicted in 
FigH are neglected to a great extent. It is natural to 
expect that such fluctuations smear out the contribution 
to the spectral function from excitations which may be 
classifled as incoherent background and only the domi- 
nating quasiparticle contributions, which may be seen as 
bright patches in FiglHl are preserved in the real system. 

Calculations based on the phenomenological approach 
to disordered charge stripes and antiphase spin domains 
give rise to a pattern formed by regions of high spectral 
intensity in the IBZ which strongly resembles ARPES 
spectra [24. .281] . Unfortunately no microscopic justifica- 
tion has been provided, of phenomenological one-body 
Hamiltonians which has been applied to derive the spec- 
tral density by means of calculations based on this sce- 
nario. Our calculation is based on the microscopic t-t'- 
t"-J model. On the other hand it seems that disorder 
may give rise spreading of spectral weight over the whole 
anti-nodal region. Such a spreading is not observed in the 
results of our calculation. In a previous paper we have 
demonstrated that the magnetic structure of the stripe 
which we have considered here, is likely to have lowest 
energy at and above the doping level 12.5%, if the dis- 
tance between axes of nearest stripes is 4 lattice spacings, 
as it has been suggested by experiments [s^ ■ Bond order 
parallel to stripe axes and long range AF order coexist in 
this magnetic structure. It has been suggested that the 
shape of the Fermi surface seen in ARPES spectra from 
Nd-LSCO and LSCO at the doping level in the range 
12.5% — 15% may be also explained in the framework of 
more conventional band calculations which neglect the 
formation of nano-scale inhomogeneities It seems 

to be hard to reconcile such a way of thinking with the 
evidence for stripes forming in these systems. 

In conclusion, motivated by results of a previous cal- 
culation indicating that, at the doping level 1/8 and 
above, the stripe structure, which consists of a) hole- 
filled two-leg ladder-like DWs with the spin-Peierls order 
formed on legs and b) AF domains of width 2 lattice 
spacings and with the changing phase of the sublattice 
magnetization by tt across each DW, is stable, we have 
performed the calculation of the single-particle spectral 
density which is generated in such a system. Our analysis 



has been made in the framework of the t-t'-t"-J model, 
with parameter values in the range suggested by compar- 
ison between band structure calculations and the Fermi 
surface of overdoped LSCO systems. The calculation 
which we have performed is a combination of the bond 
fermion method and the spin polaron approach. We ob- 
serve pronounced spectral weight both in the anti-nodal 
and nodal regions. Very similar features may be seen in 
ARPES spectra from LSCO and Nd-LSCO at the filling 
level 15%, which is exactly the same as we have assumed 
in the calculation. This similarity is not trivial, because 
different optional structures of bond order in DWs, as 
the spin-Pcicrls order on rungs in the ladder-like DWs 
or the bond order which takes the shape of two layers 
in a brick-wall, give rise to spectra at the Fermi level 
which have completely different forms. We consider the 
observed agreement between experiment and theory as 
an argument for the scenario of coexisting bond and long 
range AF orders in the stripe phase of doped cuprates. 
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APPENDIX 

We list here contributions to matrix elements H39|l . 
The contribution to (|39|l is uq for the following values 
of indices (z',/), (nj), (2,1), (0,0), (2,1); (3,0), 

(0,0), (3,0); (4,0), (0,0), (4,0); (5,1), (0,0), (5,1). The 
contribution is for (0,0), (0,0), (0,0); (0,1), (0,0), 

(0,1); (1,0), (0,0), (1,0); (1,1), (0,0), (1,1); (4,0), (0,0), 

(4.0) ; (4,1), (0,0), (4,1); (5,0), (0,0), (5,0); (5,1), (0,0), 

(5.1) ; for (0,0), (0,0), (0,0); (1,0), (0,0), (1,0); (4,1), 

(0,0), (4,1); (5,1), (0,0), (5,1); ^ for (0,1), (0,0), (0,1); 

(1,1), (0,0), (1,1); (4,0), (0,0), (4,0); (5,0), (0,0), (5,0); 

for (2,0), (0,0), (2,1); (2,0), (0,0), (3,0); (2,0), (0,-1), 
(2,1); (3,1), (0,0), (2,1); (3,1), (0,0), (3,0); (3,1), (0,1), 
(3,0); (6,1), (0,0), (6,0); (6,1), (0,0), (7,1); (6,1), (0,1), 

(6.0) ; (7,0), (0,0), (6,0); (7,0), (0,0), (7,1); (7,0), (0,-1), 

(7.1) ; for (2,1), (0,0), (3,0); (2,1), (0,1), (3,0); (3,0), 
(0,0), (2,1); (3,0), (0,-1), (2,1); (6,0), (0,0), (7,1); (6,0), 
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(7.1) , (0,1), (7,1); (7,1), (0,-1), (7,1). 
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FIG. 1: Elementary cell of the underlying spin structure as- 
sumed in the calculation (inside the dashed rectangle). Ovals 
represent singlets. 
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FIG. 2: The mechanism which gives rise to FM coupling be- 
tween nearest spins which belong to different domains. No- 
tice, that these spins would belong to the same AF sublattice 
if the system were homogeneously ordered. 
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FIG. 3: Graphical representation of some states involved in 
the process of hole propagation inside an AF domain. Zig- 
zag lines represent "broken bonds" contribution from which 
to Ising energy is higher by J/2 than contribution from bonds 
occupied by two antiparallel spins. 
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FIG. 4: Illustration of processes contributing to some terms 
in the effective Hamiltonian. These terms define the on-site 
energy and the quasiparticle hopping between sites which be- 
long to AF domains. 
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FIG. 5: Graphical representation of states involved in gen- 
eration of quantum fluctuations in |n). Contributions from 
these fluctuations determine to a large extent the shape of the 
spectral function. 
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FIG. 6: Band structure in the stripe phase, obtained within 
the scenario of coexisting bond and AF orders. 
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FIG. 7: Spectral weight intensity A~(k,Lj) along some direc- 
tions at the doping level 15%. Lorentzian broadening with 
the width O.lt has been applied. 




FIG. 8: Intensity map of the spectral weight A (k,a;) at Ef- 



